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ABSTRACT 

Polarized far-infrared and submillimeter emission is calculated for models of nonspherical dust grains 
that are constrained to reproduce the observed wavelength-dependent extinction and polarization of 
starlight. For emission from regions where the magnetic field is perpendicular to the line-of-sight, the 
far-infrared emission is expected to have substantial linear polarization at wavelengths A > 100 ^m, 
but the degree of linear polarization, and its variation with wavelength, is model-dependent. Models 
in which the starlight polarization is produced by both amorphous silicate and graphite grains have 
linear polarizations between 6% and 10% at A > 100 /xm, but for some models in which only silicate 
grains are spheroidal, the linear polarization increases from about 3% at 100 /xm to about 15% at 
1 mm. We briefly discuss the implications of these results for removal of the polarized dust emission 
from maps of the Cosmic Microwave Background, as well as the possibility of discriminating among 
interstellar dust models based on observations of far-infrared and submillimeter linear polarization. 
Subject headings: radiation mechanisms: thermal — dust, extinction — infrared: ISM 
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1. INTRODUCTION 

The presence of dust grains in the interstellar medium 
has been apparent since the recognition that obscuring 
materi al existed in at least some regions of interstellar 
space (|Barnardlll907l llQlOt lTrumpleilll930[ ). Our knowl- 
edge of interstellar dust has advanced considerably since 
then, b ut remains incomplete [see, e.g . , the review by 
iDraing f2003f ). or the book by [WhitteJ (|2Q0l)] . 

In the Milky Way, large fractions of abundant elements, 
such as Si and Fe, appear to have been depleted from 
the gas phase ( ) Jenkins! l2004l ). with the missing atoms 
believed to be present in dust grains. Studies of interstel- 
lar "reddening" produced by wavelength-dependent ex- 
tinction revealed broad spectroscopic features at 9.8 /im 
and 18 ^m, attributed to amorphous silicate material in 
submicron particles, and a strong extinction "bump" at 
2175 A, usually attributed to aromatic carbon, perhaps 
graphite or polycyclic aromatic hydrocarbons (PAHs), in 
particles with radii a < 0.01 /im. Another broad extinc- 
tion feature, at 3.4 /um, is typical of C-H stretching modes 
in hydrocarbons. Finally, the observed infrared emission 
from dust includes spectral features at 3.3, 6.2, 7.7, 8.6, 
11.3 and 12.7 /im that can be explained by a substantial 
population of particles with PAH composition and sizes 
extending down to about 40 atoms. 

These spectroscopic features indicate that amorphous 
silicates and carbonaceous materials, including PAHs, 
contribute a substantial fraction of the interstellar grain 
mass, and models developed to try to reproduce the ob- 
served extinction using observed abundances of the el- 
ements invariably employ amorphous silicates and car- 
bonaceous materials as their ma in components. 

The serendipitous discovery (Hiltner 1949': 'Halll |194l ) 
of the linear polarization of starlight demonstrated that 
interstellar grains are both appreciably nonspherical and 
aligned with the interstellar magnetic field, so that the 
dusty interstellar medium acts as a polarizer. Studies 
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of starlight polarization in the ultraviolet (IClavton et al.l 
Il992t [Anderson et al.lll996l : iMartin et al.lll999| ) indicate 
that the smallest dust grains are apparently minimally 
aligned (Kim & Martin 1995), and studies in the infrared 
[see, e.g.. IWhittetl (|2003i) ] find strong polarization in the 
9.8 /im Si-0 stretch absorption feature, requiring that a 
substantial fraction of the amorphous silicate particles be 
aligned. Recent measurements of the polarization in the 
neig hborhood of the 3. 4 /im C-H stretch absorption fea- 
ture dCh lar "erani200i l find only upper limits on possible 
polarization associated with the feature, indicating that 
the hydrocarbon material responsible for this feature is 
primarily in a grain population that is not appreciably 
aligned, which may rule out models where the hydrocar- 
bon material is located in mantles coating silicate cores. 

While we know for sure that interstellar grains are not 
spherically symmetric, their actual shapes are uncertain. 
Because the degree of alignment is unknown, we can only 
say that the grains must depart from spherical symmetry 
sufficiently to produce the observed polarization. 

Interstellar grains radiate strongly in the far-infrared. 
Because the grains are nonspherical and partially alig- 
ned, this emission must be partially polarized. The 
polarization of the far-infrared emission has been mea- 
sured for infrared-bright regions in star-forming molec- 
ular clouds (see the review by Hlldebrand 2005). The 
degree of polarization of the thermal emission from 
dust in these regions is observed to have a compli- 
cated wavelength dependence, with the polarization at 
350 /im observed to be significantly smaller than the 
100 /im and 850 /xm polarization ( Vaillanco urt ,2003: 
Hildebrand fc Kirbvll200l : lHildebrand 2005; V aillancourt! 
2007f). which has been interpreted as indicating two sep- 
arate populations of emitting grains in these clouds. 

The polarization of the large-scale submillimeter emis- 
sion from dust in the interstellar medium of the 
Milky Way has now been measured by s everal exper- 
iments, including the Archeops balloon (jBenoit et al.l 
!2004l: IP onthieu et al. 2005) and the WMAP satellite 
(|Page et all 120071 : IGold et aLl!2008f ). A quantitative un- 
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derstanding of this emission is important for removal of 
polarized Galactic "foregrounds" from maps of the cos- 
mic microwave background (CMB), but also offers an 
opportunity to test models for interstellar grains. 

In this paper, we construct models of interstellar 
dust that are consistent with the observed wavelength- 
dependent interstellar extinction, and the observed linear 
dichroism of the interstellar medium. We will calculate 
the polarized infrared and submillimeter emission from 
different models of interstellar dust heated by the general 
interstellar starlight, and, for each of them, compute the 
corresponding expected percentage polarization. This 
will allow us to gain a sense of the degree to which the 
predicted far- infrared polarization is model-dependent. 

We discuss the observed wavelength-dependent extinc- 
tion and polarization of starlight in fJ21 and the basics of 
the dust models we consider in ^ gl provides a descrip- 
tion of the technique we use to compute optical cross sec- 
tions, while [JS] introduces the formalism associated with 
"picket-fence" alignment. In SJHl we present the size dis- 
tributions and alignment functions that provide the best 
fit to the observed extinction and polarization discussed 
in We describe the formalism we use to calculate 
the partially-polarized far-infrared emission by interstel- 
lar dust in !j7l The corresponding results are presented 
and discussed in !J51 Finally, offers a summary of the 
main results of this paper, and a brief discussion of their 
consequences for studies of the CMB. 

2. EXTINCTION AND POLARIZATION OF STARLIGHT 

The grain models considered here will be required to 
reproduce the observed extinction by interstellar dust. 
Gas and dust are generally well-mixed, and the mixture 
is characterized by the observed extinction cross section 
at wavelength A per H nucleon 

A A /mag 



(obs).>N _ 



(A) EE 0.4 (In 10) 



(1) 



where A\ is the measured extinction, and Nn the col- 
umn density of H nucleons. In diffuse sightlines, JVh can 
be determined from ultraviolet observations of Lyman- 
a absorption by atomic H and Lyman- and Werner- 
band absorption from H2. The extinction produced 
by dust in the diffuse interstellar medium has been 
studied from about 0.1 fxni to sli g htly over 10 fiui 
[see, e.g., the review bv iFitzpatrickl (|2004D ]. At wave- 
lengths 4 < A 8 Mm , the extinction is somewhat 
controversial ([Lutz et all [19 96; Rosenthal ct al. 2000; 
llndebetouw et all l2005nFlahertv et all I2007i ). but the 
present study will be limited to the wavelength range 
0.1 < A < 2.65 /im, where the extinction by dust 
in diffuse regions appears to be securely known. We 
take the wavelength de pendence of t he ob served extinc- 
tion to be given by the IFitzpatrickl ()1999 l parametriza- 
tion for Ry =_Ay / {Ab — Ay) = 3.1, and, following 
iBohlin et all (|1978[ ). the value of 7Vh/(^b - Ay) is set 
to 5.8 X 10^^ cm~^/mag. 

The light emitted by stars is usually unpolarized, and 
the observed fractional polarization is produced by linear 
dichroism of the interstellar medium - a difference in the 
attenuation between two orthogonal linear polarization 
modes. We define the fractional polarization 



Fi — F2 
F1+F2 



tanh(AT), (2) 



where At = {t2 — ti)/2, Fi and F2 are the fluxes in each 
of two orthogonal polarization modes, and ri and T2 the 
extinctions for these modes. The observed polarization 
fraction p(A) tends to peak at a wavelength Amax between 
about 0.5 and 0.9 /im, with the larger values of Amax com- 
ing from sightlines through dense cloud s; diffuse clou d 
sighthnes tend to have Amax ~ 0.55 /Ltm (|Whittetj|2003l ). 

The linear dichroism is characterized by a polarization 
cross section at wavelength A per H nucleon 



(obs) 
''pol 



(A) 



At 
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arctanh[p(A)], 



which, for p ^ 1, can be approximated as 
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The strength of the polarization is characterized by the 
polarization per unit extinction Pmax/^(Amax)j where 
Pmax = p(Amax). On somc sightlincs the polarization per 
unit extinction is as large as Pmax/^(Amax) ~ 0.03/mag 
((Whitte t 2003, ) . Smaller values of PmaxM(Amax) are pre- 
sumed to arise where the local magnetic field is not or- 
thogonal to the sightline, where the magnetic field di- 
rection varies along the line-of-sight, or where the pro- 
cesses responsible for grain alignment are for some rea- 
son less effective. A successful grain model must be 
able to account for the largest observed values, i.e., 

PmaxM(Amax) « 0.03/mag. 

Various approaches have been taken to describ- 
ing the wavelengt h dependence of in terstellar po- 
larization [see, e.g.lMartin et al.l (|1999f )]. Following 
iDraine fc Allaf Akljaril ( 2006| ), we take it to be given by 
the empirical "Serkowski law" (jSerkowskilll973f ) 



p(A) « Pn 



exp 



-K In 



A 



An 



(5) 



for 0.15 < A < 1.39 /im, with Amax = 0.55 and 
K ~ 0.92. For 1.39 < A < 5 /im, we assume the power- 
law dependenc e p(A) oc A~^'^, consistent with the obser- 
vations bv .Martin et al.l ()1992t ). 

3. DUST MODELS 



As in previous stud i es (IWeingartner fc Draind 120011 : 
IDraine fc Allaf-Akbaril Ho^T 

we model the dust as 
a mixture of carbonaceous particles (including PAHs) 
and amorphous silicate grains. Models of the in- 
frared emission fro m interstellar dust (jLi fc Draindl200]] : 
IDraine fc Lil |2007| ) require a population of very small 
PAH particles, containing C/H w 55 ppm, or about 22% 
of all of the C atoms if the C abundance in the ISM is 
equal to the the current best estim ate for the solar abun- 
dance (C/H)q = (247 ± 28) ppm (|Asplund et a ll 12005) . 
We assume t he PAHs to have th e absorption cross sec- 
tions given bv IDraine fc Lil (|2007D . with a strong absorp- 
tion feature at 2175 A, and to contribute (Text (PAH, A) to 
the overall extinction. T he PAH population is exp ected 
to be randomly-oriented ([Lazarian fc DraiTiill999D . and 
therefore is assumed to contribute no polarization. This 
is consistent with the observed weakness or absence 
of polarization in the aforementioned 2175 A feature 
(|Clavton et al.l[T99l IWolff et al.l[T99l . 

In addition to the PAHs, the model includes pop- 
ulations of larger silicate particles and carbonaceous 
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grains. The size distribution of these particles is a pri- 
ori unknown, as is the degree of alignment. To de- 
termine these quantities, we will f ollow the approach 
pioneered by iKim fc Martini (|1995f ). seeking to find a 
model grain size distribution and size-dependent degree 
of alignment that can simultaneously reproduce both the 
wavelength-dependent extinction cross section (Text (A) 
and the wavelength-dependent polarization crpoi(A). 

In our models, the dust grains will be assumed to be 
oblate spheroids, with semimajor axis a along the sym- 
metry axis a, and b > a perpendicular to this axis. For 
a given axial ratio b/a, we characterize the grain size by 
the effective radius flofr = (a6^)^/'^, which corresponds to 
the radius of an equal-volume sphere. 

4. OPTICAL CROSS SECTIONS 

For a nonspherical grain of solid volume V, we define 
the extinction efficiency factor 



(6) 



-*off 



where Coxt is the extinction cross section, and a^s the 
radius of an equal- volume sphere, i.e., (3y/47r)^/'^. The 
extinction cross sections for polarized light are calcu- 
lated for oblate spheroids using the "extende d boundary 
condit ion method" (EBC M) introduced by IWatermaiil 
(119711) and develope d bv iMishchenko fc Travii (|1994f ) 
and IWielaard eTaLl (fl99l . We use the open-sourc e 
code aunpld . Ip . f made available bv lMishchenkol (|2000r ). 
EBCM codes encounter numerical difficulties when the 
target becomes large compared to the wavelength. For 
the silicate oblate spheroids {b/a > 1), we find that 
ampld . Ip . f sometimes fails for b/X > 3.885. Thus, for 
b/X > 3.885, and for both b/a = 1.4 and b/a = 1.6 (val- 
ues we will focus on from now on), we take 



Qc.t[b/X, b/a, X] « Qcxt [3.885, b/a, X] 



(7) 



for silicates. The EBCM assumes an isotropic dielec- 
tric function, but graphite is a highly anisotropic ma- 
terial. The eigenvalues of the graphite dielectric ten- 
sor are (e|| , e_L, ej,). In order to estimate Qext and Qpoi 
for spheroids composed of graphite, we use the so-called 
"1/3-2/3 approximation": we assume that 1/3 of the 
spheroids have an isotropic dielectric function e = e||, 
while 2/3 of them have e ~ e±. This gives the correct 
cross sections in the "elect ric dipole limit" aes <^ X, and 
iDraine fc Malhotral (|1993D used the discrete dipole ap- 
proximation to show that the 1/3-2/3 approximation is 
reasonably accurate even for finite grain size. The ac- 
tual composition of the carbonaceous grains is unknown. 
Based on the profile of the 3.4 /im C-H stretching mode, 
iPendleton fc Allamandolal ()2002f ) estimate that about 
85% of the carbon is in aromatic structures. We do not 
expect interstellar grains to be single-crystal graphite, 
but it is hoped that the approach adopted here will pro- 
vide cross sections that will approximate the carbona- 
ceous material present in the interstellar medium. 

The dielectric function for graphite is somewhat more 
extreme in the ultraviolet than that for silicate. For 
the graphite oblate spheroids, we find that ampld. Ip.f 
sometimes fails for b/X > 2.023. Thus, for graphite 
grains with b/X> 2.023, we take 



The top panel of Figure [T] shows Qcxt as a function 
of 2'Kb/X at four selected wavelengths for randomly- 
oriented^ silicate spheroids with b/a = 1.4 and b/a = 1.6. 
In the limit b/a —> oo, we expect QextTra^g to ap- 
proach twice the geometric cross section, with 50% of 
the total extinction produced by rays that impinge on 
the grain, and 50% contributed by small angle scatter- 
ing resulting from diffractio n [see, e.g., the discu ssion of 
the "extinction paradox" in iBohren fc HuffmanI (jT983)]. 
If (1/3,1/3,1/3) of the spheroids have symmetry axis 
(a II X, a II y, a || z), the extinction efficiency factor for 
light propagating in the x, y, or z directions is such that 



tib/X 



oo I 



[b/a 



,2/3 



2{a/b) 



1/3 



(9) 



?oxt[&/A, b/a, X] « Qcxt [2.023, b/a. A]. 



(8) 



The horizontal dashed line shows this limiting value; we 
see that the calculated cross sections appear to be con- 
verging to this limit. Our assumption summarized by 
equation ([7]) will overestimate Qcxt slightly. However, for 
realistic size distributions, the extinction at A > 0.1 /im 
is generally dominated by grains with b/X < 1. 

The bottom panel of Figure [1] is the same as the top 
panel, but for graphite spheroids. Note that while Qcxt 
for the silicate spheroids reaches values as large as 4.7 
(for A — 0.215 A in the top panel), the carbon grains 
reach only Qcxt ~ 3.1 for b/a = 1.4, and Qcxt ~ 2.8 
for b/a — 1.6; the refractive index matters! We also 
see that the extinction approaches twice the geometric 
cross section for b/X > 2, so we may use equation ([8|) to 
estimate Qcxt for large values of 6/ A. 

Finally, Figure [2] shows extinction cross sections as a 
function of wavelength for silicate and carbon spheroids 
of selected values of the effective radius Ocff. 

5. PICKET-FENCE ALIGNMENT 

The pr ocess of interstellar grain alignment [see, e.g., 
iRobergd (|2004)] leads to the short axis of grains (the 
axis with the largest moment of inertia) tending to align 
with the magnetic field direction. For oblate spheroids, 
the symmetry axis a therefore aligns with the magnetic 
field, and maximum polarization of starlight occurs for 
sightlines perpendicular to the magnetic field direction. 

We adopt a coordinate system with the line-of-sight in 
the z direction, and the magnetic field in the x direction. 
For simplicity, we will assume "picket-fence" alignment, 
i.e., that a fraction (1 -I- 2/) /3 of the oblate spheroids are 
oriented with the symmetry axis a in the x direction, a 
fraction (1 — /)/3 with a || and a fraction (1 — /)/3 with 
a II z. Thus, / = corresponds to random orientation 
(1/3 of the grains in each orientation), and / = 1 to 
perfect alignment (all of the grains with a || x). 

"Picket-fence" alignment is not physically realistic. In- 
terstellar grains will have a continuous distribution of 
orientations, and a grain with angular momentum that 
is not aligned with B will precess around B. The advan- 
tage of assuming picket-fence alignment is that for each 
grain size and wavelength, the extinction cross section 
needs to be calculated for only 3 grain orientations. We 
expect that picket-fence alignment with alignment frac- 
tion / will closely reproduce the polarization produced 

^ Random orientation is approximated by averaging over the 
three independent orientations of spheroids in the "picket-fence" 
approximation described in JS] 
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Fig. 1. — Qoxt(fe/A, b/a, A) vs. 27rf)/A at four wavelengths for randomly-oriented silicate (top) and graphite {bottom) spheroids with axial 
ratios b/a of 1.4 and 1.6. The vertical dashed lines show tiic maximum value of 2iTb/\ allowed by the EBCM calculations with the silicate 
dielectric function. Horizontal dashed lines indicate Qext{b/X oo,6/a, A) in each case. 



by grains with a continuous distribution of alignments 
provided /= (3/2) [((a • x)^) - (1/3)] . 

Let Ce:x.t,j{'m,(ieS,i>/a, X) be the extinction cross sec- 
tion at wavelength A for a spheroid of material m, with 
effective radius and axial ratio b/a, where j refers to 
the orientation of the grain relative to incident polarized 
light. We discrctizc the size range into N size bins aeff,i 
(with i e [1,7V]), uniformly spaced in log(aeff,i), with 
nsii,i and ricar,* the numbers of silicate and carbonaceous 
grains per H nucleon, and /j the fractional alignment 



of grains of size aes,i- The process of interstellar grain 
alignment remains uncertain, and it is possible that the 
alignment function / might depend on grain composition 
as well as on grain size. In this paper, we will assume 
that fi depends only on grain size, not on the material. 

The "average" extinction curve we will fit to is an av- 
erage over sightlines making various angles to the mag- 
netic field direction, so it is appropriate to compare the 
observed extinction with the extinction calculated for 
randomly-oriented spheroids. 



I 

We define the dimensionless efficiency factor for extinction by randomly-oriented grains through 
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Fig. 2. — Extinction cross sections Qcxt (leff i V^^i /'^cff vs. inverse wavelength 1/A for (leji) silicate, and (right) graphite spheres and 
spheroids with axial ratios of b/a = 1.4 and b/a = 1.6. Curves are labeled by a^Q (in /im). 
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Fig. 3. — AQpoi/aeff vs. inverse wavelength for [left) silicate spheroids, and (right) graphite spheroids with axial ratios b/a = 1.4 and 
b/a = 1.6. Curves are labeled by a^ff (in /im). The pronounced drop in Qpd for graphite grains as A increases from 0.20 /im to 0.22 /ivci is 
a result of the wavelength dependence of the dielectric function ex (see text). 

The index j runs from 1 to 3 as we require three orientations for picket-fence alignment with oblate spheroids: 

j = 1 for k ± a and E || a, 

j = 2 for k ± a and E ± a, (11) 
j = 3 for k II a and E ± a, 

where k is the incident wavevector. We also define an efficiency factor for polarization by perfectly-aligned grains: 

aeff 6\ 11 



Qpoi m, 



A a 



(12) 



Figure [3] shows the polarization efficiency factors Qpoi as a function of wavelength A for selected sizes of silicate and 
carbon spheroids. Note that Qpoi tends to peak near 27raoff/A « 1.6 for the silicate spheroids, whereas it peaks near 0.8 
for the carbon spheroids. Small carbon spheroids show a strong peak in Qpoi near A « 0.23 /im; this is a result of the 
strong peak in Im(ex) due to tt — > tt* electronic transitions near this frequency. 

The model we consider here involves a mixture of silicate and carbon (including PAHs) dust grains. In such a case, 
the extinction cross section per H nucleon is given by 



N 



(A) (PAH,A)+^^ 



m, ■ 



(13) 
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where we sum over materials m = {sil, car} and over grain size a(.s,i, with nm,i the number of grains of composition 
m and size ace,i per H nuclcon. The polarization cross section for the model is 

N 

<Jpo\{X) — ^ ^ fi nrn,i Tra^ff J Qpoi 

m i—1 

A priori, we do not know the grain shape, except that the grains must be sufficiently nonspherical to be able to 
reproduce the observed polarization of starlight. We will proceed using trial grain shapes, and will see whether we are 
able to reproduce the polarization with a size-dependent alignment fraction /(a) < 1. 

For a given model (with assumed axial ratio (6/a)„i for material m), we use the technique described in |2]to obtain the 
cross sections Cextj [w, aeff,i, (6/a)m, A], and then seek the abundances n^^m and alignment function fi that minimize 
a positive-definite penalty function 

^ = E(*^)'' (15) 

t 

where the ^f 's are measures of deviation from an ideal fit to the data. We use the penalty function described by 
iDraine fc Allaf-Akbaril (|2006D . which includes a contribution proportional to the square of the deviation of the model 
extinction curve from the observed extinction curve, and a contribution proportional to the square of the deviation of 
the model polarization curve from the observed polarization curve, as well as terms that measure the non-smoothness 
of the size distributions Usn^i and ncar.i, and of the alignment function fi. The degree to which the amount of material 
in the model grains deviate s from the am ount of material "depleted" from the gas, as inferred from observations of 
ultraviolet absorption lines (jJenkinsll2004l ). also contributes to the penalty P. 



m, 



acS-, 



(14) 



TABLE 1 

Models 



Model 


b/a 


b/a 


Vsil 


V ^ 

V car 


number 


(silicate) 


(carbon) 


(10-27cmVH) 


(10-27cm3/H) 





1.0 


1.0 


3.00 


1.91 


1 


1.4 


1.0 


3.81 


1.93 


2 


1.4 


1.4 


3.34 


1.89 


3 


1.6 


1.0 


3.18 


1.92 


4 


1.6 


1.6 


3.20 


1.89 



"Observed" 2.50 1.00 



''Including C/H = 55 ppm in PAH material, with assumed carbon density 2.24 g/cm^. 



6. SIZE DISTRIBUTIONS AND ALIGNMENT FUNCTIONS 

Polariz ation in th e 9.7 /im silicate feature is routinely 
observed (jSmith et al. 2000), and therefore it is clear that 
the silicate grains must be nonspherical. As a result, all 
of our models (1-4) include spheroidal silicate grains, 
with axial ratios b/a — 1.4 or b/a — 1.6. Because there 
is no observational evidence of polarization produced by 
carbonaceous grains, we will consider two possibilities, 
namely, that the carbonaceous grains are spherical, or 
just as nonspherical (and ahgned) as the silicates. 

We study four models involving nonspherical grains. 
Model 1 has spheroidal silicate grains with b/a = 1.4, 
and spherical carbon grains. Model 2 assumes both sil- 
icate and carbon grains to be spheroids with b/a — 1.4. 
Model 3 considers silicate spheroids with a more extreme 
axial ratio, b/a — 1.6, but with spherical carbon grains. 
Finally, our fourth model assumes both silicate and car- 
bon grains to be spheroids with b/a = 1.6. 

We obtained the best-fit grain abundances n^^car and 
rij^sii, and degree of alignment /i, for these four dis- 
tinct models using the Levenberg-Marquardt algorithm 
to minimize the penalty function P discussed in ^JSl 

Figure m shows the "observed extinction" that we seek 
to reproduce, together with the extinction obtained for 
each model. For comparison purposes, we also carried 



out a fit to the extinction using spherical silicate and car- 
bon grains (model 0), with no consideration of polariza- 
tion. Note that all models are more-or-less equally good 
at reproducing the observed extinction. All of the mod- 
els tend to underpredict the extinction for A > 1.2 /im, 
suggesting that perhaps our adopted composition is in- 
sufficiently absorptive at A 2 /um. It should be kept in 
mind, however, that the "observed" extinction is some- 
what uncertain at the longer wavelengths, owing to the 
increasing difficulty in determining the extinction per 
H nucleon when the H column density cannot be de- 
termined direcU^Y_jasmgLynm absorption measure- 
ments. 2 iNishivama et all (|2006f) find Ax cx A-^-^^ for 
1.25 < A < 2.2 /im, steeper than the "observed" extinc- 
tion used here, and closer to the behavior in our models. 
The volumes of silicate and carbonaceous materials re- 
quired by each model is given in Table [1] 

Elemental Abundances 

The most recent estima tes of solar ab undances indicate 
that Mg/H = 38 ppm (jGrevesse k Sauval .1998.). an d 
Si/H = 32 ppm and Fe/H = 29 ppm (|As2lun3l2OOO0 . 



^ To accurately measure the weaker extinction in the infrared, 
larger dust columns are needed, but the ultraviolet extinction is 
then so strong that Lyman-a measurements become infeasiblc. 
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Fig. 5. — Mass distribution for (ie/t) silicate and (right) graphite spheres a nd sp heroids. Results are shown for models through 4. The 
curves labeled "DL07" show the size distributions of iWeingartner fc Drain3 l|2001l) . as modified bv iDraine fc: Lii (i2007i) . 



If all of the Si is consumed, we could form silicate ma- 
terial with the olivine-like composition Mgi.iFeo. 98104. 
The density of Si atoms in fayalite (Fe2Si04) is about 
1.3 X 10^2 cm-3. With this density of Si, Mg^ ;iFeo.9Si04 
would have density p « 3.7 g/cm.^, and the solar abun- 
dance of Si could produce V = 2.5 x 10~^^cm'^/II of 
silicate material. Thus, we see that the silicate volumes 
in Table [T] exceed what appears to be allowed by current 
estimated solar abundances by 20% to 50%. 

For a carbon density of 2.2 g/cm^, the carbon volume 
in model corresponds to C/H = 210 ppm. On the 
well-studied sightline toward C Op h, the gas-phase car - 
bon abundance is about 135 ppm (iCardein et al.lll993l) . 
Adding this to our estimate for the carbon in dust grains 
would give an interstellar C/H = 345 ppm, 40% greater 
than the most recent est imate of the solar a bundance 
(C/H)o = 245 ± 30 ppm (lAsplund et al.l[200l . 

The amount of grain material called for in these fits is 



cause for concern - our models tend to require about 40% 
more silicates and carbon than appears to be available if 
interstellar abundances are the same as solar abundances. 
However, interstellar abundances of C, Mg, Si, and Fe 
may exceed the abundances of these elements in the Sun. 
Furthermore, solar abundances may have been underes- 
timated. For example, it has recently been argued that 
the oxygen abundance in the solar pho tosphere has been 
underestimated by a factor of ab out 2 (jLandi et al.ll2007l : 
ICenteno fc Socas-Nava^l2008f) . 

Figure [5] shows the size distributions for the silicate 
grains and the carbon grains for our five dust models. 
Two classes of models are clearly visible. The models 
where only the silicate grains are aligned (models 1 and 3) 
have broad distributions for the silicates, while the car- 
bonaceous grains have a pronounced peak near 0.08 /im. 
On the other hand, in models where both silicate and 
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carbonaceous grains are aligned (models 2 and 4), the 
silicate size distributions have a dominant peak near 
0.15 ^m, while the carbonaceous grains have pronounced 
peaks near 0.03 ^m and 0.25 /xm. Despite these differ- 
ences, all of the models 1 through 4 produce polarization 
that is in reasonable agreement with the observations, as 
can be seen in Figure [6l but this is achieved with quite 
different alignment functions /(a), also shown in Fig- 
ure [SI The models with spherical carbon grains (mod- 
els 1 and 3) have alignment functions that rise rapidly 
from near-zero to near-unity at a « 0.1 fim. When non- 



spherical carbon grains are included, however, the large 
polarization cross sections for the carbon grains mean 
that the required polarization can be achieved without 
complete alignment of the a > 0.15 /um grains. For 
model 2 {b/a — 1.4), the alignment for a^B ^ 0.1 /im 
is / w 0.4; when 6/a is increased to 1.6, the alignment 
of the larger grains drops to about 0.3. Note that these 
alignment functions reproduce the highest observed val- 
ues of Pmax/A-v- It is possible that, in some regions, the 
dust grains have lower degrees of alignment. 



7. FAR-INFRARED CROSS SECTIONS AND EMISSIVITIES 

At wavelengths A ^ a (the "Rayleigh limit" ) , the absorption cross section of a grain depends only on the orientation 
of the oscillating electric field relative to the grain shape, and not on the direction of the oscillating magnetic field 
(i.e.. Cabs, 3 = Cabs, 2)- For a spheroid, the absorption cross section is then 



Cabs = Cabs.l cos C + Cabs,2 sin (, 



(16) 



where ( is the angle between the electric field E and the symmetry axis a. Randomly-oriented spheroids have 
{cos^ () = 1/3, which leads to the "average" absorption cross section 



(a 



abs/ 



2 (Cabs, 1 



+ 2 Cabs, 2)- 



(17) 



Let 0Q be the angle between the direction of the magnetic field, x, and the direction of propagation. We denote quan- 
tities for radiation with E parallel (respectively, perpendicular) to the (k, x) plane with the label || (respectively, _L) 
If a fraction / of the grains is perfectly aligned, and a fraction (1 — /) is randomly-oriented, 

1-/ 



cos 



C||) = 

cos'-C±) = 



3 

1-/ 
3 ■ 



/sin^ 



The polarization-averaged absorption cross section is then given by 



c 



abs 



2 (^abs 



c'l ) 



a 



abs,l 



§-/(cos2' 



c 



abs, 2 



o + / cos' 



^°-3 



whereas the cross section for polarized emission is defined as 



c 



1 



pol 



abs 



c 



1 



abs/ 



/ (C'abs,2 — C'abs,l) sin^ Oq. 



(18) 
(19) 

(20) 
(21) 
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Therefore, for a mixture of partially-aligned grains, the total and polarized emissivities are 

+ i (^C0S2 Oo-^^Y.j d«<=ff ^^^"^^^^^^"'^'^ - ^^''^''^ (B^)m,a.n^ (22) 

1 1 /" dn 

^ /A.poi(eo) - 2 ^0 E 7 d^-ff /(acff)(abs,2 - abs,i) (i^A)™,,,,, (23) 

(i^A).„,,../dT -^^^°-;;;.-'^) i..(T), (24) 



where 



dP(aoff, m, U)/dT is the temperature distribution for grains of type m and size OeS i n a radiation field of intensity U , 
and Bx{T) the Planck function. We use the temperature distributions computed bv lDraine fc Lil (2007) for spherical 
grains. For the spheroidal components of our models, we correct these distributions to account for the non-sphericity 
of the grains by requiring that the power they absorb equals the power they re-radiate. The correction is very small. 



8. RESULTS AND DISCUSSION 

Figure [7] shows the emission spectra calculated for each 
of the five grain models listed in Table[T]when illuminated 
by the standard interstellar radiation field for the solar 
neighborhood. Also shown is t he emission spectru m cal- 
culated for the grain model of lDraine fc H (|2007l ). The 
emission is calculated using the cross section (Cabs) 
from equation ([T7]), appropriate for randomly-oriented 
dust. These spectra of the thermal emission do not 
include rotational emission from spinning dust, which 
appears to dominate for A > 4000 um (i^ < 75 GHz) 
(jPraine fc Lazarianlll998HFinkbeineill2004D . 

The emission spectra for all six models are relatively 
similar. In each case, A (Ix) peaks near 130 /zm, with an 
amplitude within 20% of the lDraine fc Lil ()2007D (there- 
after, DL07) model. Moreover, all models agree ex- 
tremely well for wavelengths below about 20 /im, where 
the signal is entirely dominated by PAH features. 

However, there are some noticeable differences, espe- 
cially in the ratio of the 130 fim to the 24 fim emis- 
sions. Although these ratios are within 15% of each 
other for model and the DL07 model, models involv- 
ing spheroidal grains lead to smaller values of this ratio, 
typically by a factor of 1.5. This difference arises mostly 
from variations in the size distributions between models, 
with differences in grain shape also having a small effect 
on the average emission spectrum. All other differences 
are fairly modest (usually at the 10% to 20% level), and 
generally smaller than the actual uncertainties in the ap- 
propriate dielectric functions. 

The fact that all of these spectra are in good agree- 
ment overall, but with significant well localized differ- 
ences, opens a very interesting window through which 
one can hope to strongly constrain dust models with ob- 
servations of the total dust emission at just a couple of 
wavelengths. We have not attempted this exercise here. 

Figure [5] shows the polarized emission spectra calcu- 
lated for models 1 through 4, as well as the correspond- 
ing degrees of polarization (in %), for lines of sight per- 
pendicular to the direction of the local magnetic field, 
and the dust exhibiting the highest degrees of alignment 
found in the diffuse interstellar medium. 



Contrary to the total emission spectra shown on Fig- 
ure [71 the polarized emission spectra clearly reflects the 
fact that we consider two classes of models. The spectra 
of polarized emission for models 2 and 4, with spheroidal 
silicate and carbonaceous grains, are indeed remarkably 
similar, but quite different from those for models 1 and 3, 
with silicate spheroids and carbonaceous spheres. This 
translates into interestingly different predictions for the 
expected degree of polarization for emission by interstel- 
lar dust, as can be seen in the right panel of Figure [5] 

Although all models predict virtually no polarization 
below 40 /im (as small grains and PAHs are not aligned) , 
the wavelength dependence of the degree of polariza- 
tion for A > 40 iim is quite model-dependent. Models 
with both silicate and carbonaceous spheroidal grains 
lead to the degree of polarization sharply increasing be- 
tween about 40 /zm and 600 ^m, after which it reaches a 
plateau. However, the degree of polarization for models 
in which only the silicate grains are spheroids increases 
steadily with increasing wavelength, without converging 
to an upper limit. As a result, models 1 and 3 predict 
significantly higher degrees of polarization than models 
2 and 4 for wavelengths longer than about 200 /im. For 
example, when A = 3000 /^m, the maximum degree of po- 
larization predicted by models 1 and 3 is of order 10%, 
while models 2 and 4 lead to a maximum of about 15%. 
These differences provide another way of discriminating 
between dust models, this time from observations of the 
polarized emission from interstellar dust. 

The strong wavelength-dependence of the polarization 
in models 1 and 3 can be easily understood. In these two 
models, the carbonaceous grains are not aligned. For 
the adopted grain optical properties, the silicate grains 
contribute an increasing fraction of the emission as the 
wavelength increases, in part because the silicate grains 
are slightly cooler than the carbonaceous grains (this ex- 
plains the behavior at A < 200 /im) , and in part because 
the ratio of the silicate opacity to the graphite opacity 
increases with increasing wavelength for A > 100 /im. 

9. SUMMARY AND FURTHER DISCUSSION 

The main results of this paper are as follows: 

1. We present models of interstellar dust based on 
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A (yLim) 

Fig. 7. — Total emission spectra for the different grain models listed in Table [T] as well as for the model described in lDraine fc Lil |2007), 
each heated by the interstellar radiation field estimated bv lMathis et al.1 l l 19831) . The x- axis is now labeled by wavelengths, an d the symbols 
correspond to the DIRBE results for HI correlated emission at high Galactic latitudes IIDw ck ct al. 1997: [A"rendt et al.lll993) . Frequencies 
covering the full range of Planck HFI frequencies (100 GHz, 217 GHz, and 857 GHz) are also indicated. 



857 217 100 GHz 657 217 100 GHz 




A (f^m) A (jum) 

Fig. 8. — (Left) Spectra of polarized thermal emission predicted by models 1 through 4 in directions perpendicular to the direction of 
the local magnetic field, and (Right) corresponding degrees of polarization for each model. Two classes of models are clearly visible. Note 
that these models do not include rotational emission from spinning dust, which will be important for A > 4000 fim (u < 75 GHz). 



mixtures of spheroidal amorphous silicate grains, 
spheroidal graphite grains, and PAH particles, with 
a degree of alignment / that is a function of the ef- 
fective radius Ocff of the grains. The models differ 
in assumptions concerning the shapes (axial ratios 
b/a = 1.4 and 1.6 are considered for the spheroidal 
grains) of the silicate and graphite grains. 

2. Size distributions dn/dacs and alignment func- 
tion /(aoff) are found that reproduce the observed 
wavelength-dependent extinction a nd polarization 
of starlight. As found previously bv lKim fc Martini 



(fl995l ). the de gree of alignment is small for 
Oeff < 0.05 /im, and large for Qcs ^0.1 /im. 

3. For models where only the silicate grains are 
aligned, the degree of alignment approaches 100% 
for a > 0.1 /im, whereas when both silicate and 
graphite grains are partially-aligned spheroids with 
b/a = 1.4, the observed starlight polarization can 
be reproduced with / ^ 0.4 for a > 0.1 /im. If the 
axial ratio is increased to b/a = 1.6, the required 
degree of alignment drops to / « 0.3. 

4. For each model, we calculate the partially-polarized 
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infrared emission when the dust is illuminated 
by the local starlight background. The total in- 
frared emission spectra are generally similar, and 
in approximate agreement with observations of the 
Galactic cirrus emission at |6| > 20°. In detail, the 
spectra are model-dependent, particularly between 
about 20 fxm and 160 /xm. No attempt has been 
made to adjust our models to reproduce observed 
FIR spectral energy distributions. 

5. Because the small grains and PAHs are not aligned, 
the polarized emission is very weak for A < 40 fim. 
However, for A > 60 /xm, the linear polariza- 
tion of the thermal emission is appreciable when 
viewed from directions perpendicular to the lo- 
cal magnetic field. At 100 fim, the different 
models considered here result in linear polariza- 
tions between about 3% and 6%. For wavelengths 
longer than 500 fim {v < 600 GHz), of relevance 
to CMB studies, the degrees of polarization pre- 
dicted by our models range from about 9% to 
around 15% for the thermal emission. However, 
for A > 4000 /im < 75 GHz), one should take 
into account the additional emission from spin- 
ning dust, which is not included here. This ro- 
tational eniissi£n_js_exgeeted_to be nearly unpo- 
larized ([Lazarian &: Drain^l2000[ ). which will lower 
the overall polarization. 

6. For models that employ more than one type of 
grain, which is the case of all our models, their 
relative contribution to the far-infrared and sub- 
millimeter emission will in general vary with wave- 
length. If the two types of grain have different 
values of Cpoi/Cabs because of different alignment 
fractions or grain properties, the degree of polariza- 
tion will vary with wavelength. This effect is very 
pronounced in our models 1 and 3 where only the 



silicate grains are aligned with the magnetic field. 

These last two results are particularly interesting for 
studies of the polarized components of the CMB. Polar- 
ized dust emission is indeed one of the brightest Galactic 
foregrounds at microwave frequencies, and being able to 
remove it from Galactic maps is one of the biggest chal- 
lenges that will limit our ability to detect B-modes. 

The maximum degrees of polarization we calculated 
in fJH] are consistent with the level inferred from th e 
WMAP observations (|Page et al.ll2007t iGold et al.ll2008f) , 
although 5% is on the lower side of our predictions. 
Upcoming experiments, such as Planck, might therefore 
want to allow for the possibility of dust polarization being 
slightly higher in their sky models, and consider how this 
would affect their ability to constrain E- and B-modes. 

More important is the possible wavelength-dependence 
of the degree of polarized dust emission. Models involv- 
ing both silicate and carbonaceous spheroidal grains lead 
to a roughly constant level of polarization at all Planck 
frequencies. However, there is currently no evidence that 
carbonaceous dust grains contribute to the polarized sig- 
nal at any wavelength, and models where only silicate 
grains are aligned spheroids lead to degrees of polariza- 
tion that are strongly frequency-dependent. Although a 
possible challenge for the CMB community, this will al- 
low the use of Planck observations to test whether the 
submm emission arises from a single kind of grain or from 
multiple grain types having different values of Cpoi/Cabs- 
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